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Abstract 

We present the first two leading terms of tlie (genus) expansion of the free 
energy for ensembles of normal and complex random matrices. The results are 
expressed through the support of eigenvalues (assumed to be a connected domain 
in the complex plane). In particular, the subleading (genus- 1) term is given by the 
regularized determinant of the Laplace operator in the complementary domain with 
the Dirichlet boundary conditions. An explicit expression of the genus expansion 
through harmonic moments of the domain gives some new representations of the 
mathematical objects related to the Dirichlet boundary problem, conformal analysis 
and spectral geometry. 
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1 Introduction 



Ensembles of random matrices have numerous important applications in physics and 
mathematics ranging from energy levels of nuclei to number theory. An important in- 
formation is encoded in the expansion {N is the size of the matrix) of different 
expectation values in the ensemble. Many relevant references can be found in [1] . 

In this paper we discuss 1/A^-expansion in statistical ensembles of normal and complex 
matrices. A matrix M is called normal if it commutes with its Hermitian conjugate: 
[M, = 0, so both matrices can be diagonahzed simultaneously. Eigenvalues of a 
normal matrix are complex numbers. The statistical weight 

of the normal matrix ensemble is specified by a potential function W{M) (which depends 
on both M and M^). Here k is a parameter, and the measure of integration over 
normal matrices is induced from the flat metric on the space of all complex matrices. 

Along the standard procedure of integration over angle variables [2], one passes to 
the joint probability distribution of eigenvalues zi, . . . ,zn. The partition function is then 
given by the integral 

1 r ^ 

Here A„(:i) = Yllj 

{zi — Zj) is the Vandermonde determinant and (Pz = dxdy for 
z = X + iy. The A^-dependent normalization factor is put here for further convenience. 

The model of normal matrices was introduced in [3]. This model is the particular 
(3=1 case of a more general one, referred to as 2D Coulomb gas with the joint probability 
distribution I A^(;2i) nf=i e^^^'^^d^Zj. 

For the potential of the form 

W{z) = -zz + V{z) + V{z) (2) 

where V{z) is an analytic function in some region of the complex plane (say, a polyno- 
mial), the normal matrix model is equivalent to the ensemble of all complex matrices 
with the same potential. It generalizes the gaussian Ginibre-Girko ensemble [4]. When 
passing to the integral over eigenvalues, the partition function for complex matrices dif- 
fers from the one for normal matrices by a normalization factor only [2]. Both models 
are then reduced to the 2D Coulomb gas (with /3 = 1) in the external potential. Note 
also a formal similarity with the model of two Hermitian matrices. Its partition function 
is given by the same formula (1), with the potential (2), but Zi and Zi are to be regarded 
as two independent real integration variables, with d'^Zi being understood as dzidzi. 

It appears that the potential of the form (2) is most important for applications [5]. 
In the main part of the paper, we concentrate on this case, so one may, in this context, 
ignore the difference between the normal and complex ensembles, taking the 2D Coulomb 
gas partition function as a starting point. Physical applications of this model include 
the quantum Hall effect, the Saffman- Taylor viscous fingering and, conjecturally, more 
general growth problems which are mathematically described as a random evolution in 
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the moduli space of complex curves. Recently, the normal matrix model was shown [6] to 
be closely related to the matrix quantum mechanics, and, therefore, to the c — 1 string 
theory. 

In addition to this it appears that the large N limit of the normal or complex random 
matrices admits a natural geometric interpretation relevant to the 2D inverse potential 
problem, the Dirichlet boundary problem and to spectral geometry of planar domains. 
In this paper we concentrate on calculation of the expansion of the free energy, 
F oc logZiv, and on its algebro- geometric meaning, leaving physical aspects for future 
publications. 

The large A'" limit also implies the hmit ^ — > 0, so that kN is finite and fixed. We 
prefer to work with the equivalent ^-expansion, rather than with the expansion, thus 
emphasizing its semiclassical nature. The free energy of the Hermitian, two-Hermitian, 
normal and complex matrix ensembles with the potential (2) has an /i-expansion of the 
form logZjv = Z]g>o where g-th term is associated with the contribution of 

diagrams with Euler characteristics 2 — 2g, in the perturbative expansion of the free 
energy. Here we discuss the first two terms, Fq and Fi. 



The leading term, Fq, is the contribution of planar diagrams, and Fi is commonly refered 
to as genus 1 correction. 

When becomes large new macroscopic structures emerge. Invoking a physical 
analogy, one may say that the gas of eigenvalues segregates into "phases" with zero and 
non-zero density separated by a very narrow interface. The domain D in the complex 
plane where the density is non-zero is called the support of eigenvalues (it may consist 
of several disconnected domains). The density at any point outside it is exponentially 
small as iV — cxD. 

The leading contribution to the free energy, the Fq term in (3), is basically the 
Coulomb energy of particles confined in the domain D. For the potential of the form 
(2) it is the tau- function of curves introduced in [7]. It encodes solutions to archetypal 
problems of complex analysis and potential theory in planar domains. 

Here we review these results and also compute the genus- 1 correction to the free 
energy. The latter is identified with the free energy of a free bosonic field in the domain 
D'^ which is complementary to the support of eigenvalues, i.e., in the domain where the 
mean density vanishes: 



Here det(— Aqc) is a properly regularized determinant of the Laplace operator in D*^ with 
Dirichlet boundary conditions. This suggests interesting links to spectral geometry of 
planar domains. 

The genus expansion in the Hermitian random matrix model beyond the leading 
order has been obtained in the seminal paper [8]. In [9], the genus-1 correction was 
interpreted in terms of bosonic field theory on a hyperelliptic Riemann surface. The 
genus 1 correction to free energy of the model of two Hermitian matrices with polynomial 
potential was found only recently [10]. 




(3) 




(4) 
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2 The planar large limit 



In this section wc briefly recall the large limit technique. This material is standard 
since early days of random matrix models (see., e.g., [11]). An appealing feature of the 
model of normal or complex matrices is a nice geometric interpretation and a direct 
relation to the inverse potential problem in two dimensions. 

As was already mentioned, the parameter h tends to zero simultaneously with N ^ oo 
in such a way that to — Nh is kept finite and fixed. Using the Coulomb gas analogy, one 
may say that the leading contribution to the free energy is equal to the extremal value 
of the energy 

E^J2^og\z,-z,\ + ^Y.W(z,) (5) 

i^j i 

Equilibrium positions of charges are given by the extremum of the plasma energy: dz^S ~ 
= 0. 

Consider the 2D Coulomb potential $(^) = —fiJ2i^og \z — z^^ created by the charges. 
Writing it as 

^{z) = - Jlog\z-C\^p{z)d^z 

where 

p{z) = -^A^z)=h^5^'\z-z,) (6) 

is the microscopic density of eigenvalues (a sum of two-dimensional delta-functions), we 
assume that $ in the limit can be treated as a continuous function. It is normalized 
as / A^{z)d'^z — —Anto. Let $o be this function for the equilibrium configuration of 
charges, then 

d,i%{z) - Wiz)) = a,($oW - Wiz)) = (7) 

with the understanding that this equation holds only for z belonging to a domain (or 
domains) where the density is nonzero. Applying to the both sides, we see that the 
equilibrium density, po{z), is equal to —^AW{z) in some domain D (the support of 
eigenvalues) and zero otherwise: 

Po(2) = I Q Ito" ^o(z) ^- J^\og\z-C\'^a(fC 

Here D*^ = C \ D is the domain complimentary to the support of eigenvalues and 

a = -j-AW{z,z) (8) 

In this and in the next section we assume the special form of the potential (2). Then 
Po = I/tt in the domain D. 

The shape of D is determined by the function V{z). Let us assume, without loss of 
generahty, that e D and parametrize V{z) by Taylor coefficients at the origin: 

V{z) = j:hz'' (9) 

k>l 
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The parameters tk (coupling constants of the matrix model) are in general complex 
numbers. Multiplying (7) by z~'' and integrating over the boundary of D, we conclude 
that the domain D is such that —TrktkS are moments of its complement, D'^, with respect 
to the functions z~^: 

tk^- \ ( z-HH =77^/ ^"'^^■^ (10) 
tt/c JD" 2mk JdD 

Besides, from the normalization condition we know that the area of D is equal to vrto- 
To find the shape of the domain from its moments and area is the subject of the inverse 
potential problem. These data determine it uniquely, at least locally. 

Here we assume that D is a connected domain. For example, in the potential W = —zz 
the eigenvalues uniformly fill the disk of radius y/HN. Small perturbations of the potential 
slightly disturb the circular shape. 

In what follows, we need some functions associated with the domain D, or rather with 
its complement, D'^. The basic one is a univalent conformal map from the exterior of the 
unit disk onto the domain D'^. Such a map exists by virtue of the Riemann mapping 
theorem. Let U be the unit disk and U'^ its complement, i.e., the exterior of the unit disk. 
Consider the conformal map z{w) from onto D*^ normalized so that z{oo) = oo and 
r = \im^^^ z{w) / w is real, then the map is unique. In general, the Laurent expansion 
of the function z{w) around infinity is 



z[w 



) =rw + J2ukW~'' (11) 



fc>0 



The real number r is called the (external) conformal radius of D. Since the map is 
conformal, all zeros and poles of the derivative z'{w) = dyjz{w) are inside the unit circle. 
We also need the function z{w) given by the Laurent series (11) with complex conjugate 
coefficients and the Green function of the Dirichlet boundary problem in D'^. In terms of 
the conformal map, the latter is given by the explicit formula 



G{z,z') ^\og 



w{z) — w{z') 



w{z)w{z') — 1 

Here w{z) is the conformal map from D*^ onto U"^ inverse to the z{w). 



(12) 



3 The leading term of the free energy 



The leading contribution to the free energy is the value of the Coulomb energy (5) 
(multiplied by h"^) for the extremal configuration of charges: 

Fq = f f \og\z - z'\a(z)a(z')dhdh' + [ W(z,z)adh 
JdJd Jd 

The integrated version of the extremum condition (7) tells us that ^o{z) —W{z) — const 
for any ^ G D. The constant can be found from the same equality at 2; = 0, and we 
obtain Fq as an explicit functional of the domain D: 



- / / log 

JDJD 



1 1 

Z z' 



a{z)a{z')dhd^z' 



(13) 
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For the special potential of the form (2), when a — 1, the free energy is to be regarded 
as a function of to and the coupling constants tk. 

Properties of Fq immediately follow from known correlation functions of the model 
in the planar large limit. See [5, 7] for normal and complex matrices and [12] for 
similar results in the context of the Hcrmitian 2-matrix model. Some of these correlation 
functions previously appeared in studies of thermal fluctuations in classical confined 
Coulomb plasma [13]. Integrable structures associated with Fq were studied in [14, 7, 15, 
16]. Here is the list of main properties of Fq for the most important case a — I/tt. 



1-st order derivatives: 



Otk TT Jd 

dt. 



(14) 



|o ^ 1 /■ 
TT Jd 



log \z\ d z 



can be combined in the generating formula 

1 



V{z)Fo = - / log \z-' - C'l'd'C , zeD', 



where 



d 



d 



d 



dto 



k>l 



dtk 



9tk, 



(15) 



(16) 



Since the derivatives of Fq with respect to the moments tk are moments of the com- 
plimentary domain, this function formally solves the 2D inverse potential problem. 



2-nd order derivatives: for z, z' e D'^ we have 

V{z)V{z')Fo = 2G{z,z') 



log 



z' 



[17) 



where G{z, z') is the Green function of the Dirichlet boundary problem in D'^ (12). 
Note that the logarithmic singularity of the Green function at z = z' cancels by 
the second term in the right hand side. In a particular case when both z, z' tend 
to infinity, we get a simple formula for the conformal radius: 



a2Fo = 21ogr 



(18) 



• 3-d order derivatives. The generating formula reads [15]: 

V{a)V{b)V{c)Fo ^-^i dnG{a,i)dnG{h,i)dnG{c,Om (19) 

An important corollary of this formula and eq. (17) is the complete symmetry of 
the expression V{a)G{b, c) with respect to all permutations of the points a, b, c. 
Another corollary of (19) is the following residue formula valid for j, k,l > [16]: 

d^Fp ^ J_ I hj{w)hkiw)hi{w) dw 
dtjdtkdti 27ri J\w\=i z'{w)z'{w~'^) w 
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Here hj{w) are polynomials in w of degree j: 

hj{w) — w-j- {z-'{'w))+ for j > 1 and hoiw) — 1 , 

(Jj UJ 

where (...)+ is the positive degree part of the Laurent series. The notation z'{w^^) 
means the derivative dz{u)/du taken at the point u = w~^. This formula is es- 
pecially useful when z'{w) is a rational function, then the integral is reduced to a 
finite sum of residues. We use this below. 

• Dispersionless Hirota equations. The function Fq obeys an infinite number of non- 
linear differential equations which are combined into the integrable hierarchy of 
dispersionless Hirota's equations. See [7, 15] for details. 

• WDVV equations. Suppose V{z) is a polynomial of m-th degree, i.e., tfc = for all 
k > m. On this subspace of parameters, Fq obeys the system of Witten-Dijkgraaf- 
Verhnde-Verlinde (WDVV) equations 

Fif-'fk = FfcF-^Fi for all < i, j, k < m - 1 (21) 

where Fj is the m by m matrix with matrix elements {?i)jk — dtot-dtk ' ^^^^ 
details. 

To conclude: Fq is a "master function" which generates objects of complex analysis 
in planar simply-connected domains. The full free energy of the matrix ensemble, F, 
may be regarded as its "quantization" . 



4 The genus 1 correction to the free energy 
4.1 The result for Fi 

We now describe the result for the genus- 1 correction Fi. We start with the special 
potential (2). Then Fi is expressed entirely in terms of the metric on the U*^ induced 
from the standard flat metric on the ;s-plane by the conformal map: dzdz = e'^'^^'^Uwdw. 
Here 

(j){w) = log \z'{w)\ (22) 

and z{w) is the conformal map U'^ — > D*^ (H)- The derivation of this formula and its 
extension to a general potential is outlined in Section 5. 

We found that 

^1 = -^/ {<l>dn(l> + 2<l>)\dw\ (23) 

247r J\w\=i 

Here 9„ is the normal derivative, with the normal vector pointing outside the unit circle. 
The derivation of this formula is outlined, for a more general model, in Section 5. 

Since (j){w) is harmonic in U'^, we may rewrite the r.h.s. of (23) as 
Fi = -^/ \V(j)\^d^w - (f (Pldwl 

247rJ|«;|>l' ^' 127rJ|«;|=l ' ' 
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Here we recognize the formula for the regularized determinant of the Laplace operator 
Aqc = 4:dzdz in D*^ with Dirichlet conditions on the boundary. The first term is the 
bulk contribution first found by Polyakov [17] (for a metric induced by a conformal map 
it reduces to a boundary integral), while the second term, computed in [18], is a net 
boundary term (see also Section 1 of [19]): 

Fi = -^logdet(-ADc) (24) 

In the particular case W{z) = —zz we get = log to that coincides with the result 
of [20] obtained by a direct calculation. 

The appearence of elements of quantum field theory in a curved space is not accidental. 
A field-theoretical derivation of this result will be given elsewhere. 



4.2 Rational case 

Before explaining the origin of the explicit formula for Fi we write it in yet another 
suggestive form. Consider a domain such that z'{w) is a rational function: 



z'{w) — r 



""-^ w - a,- 



All the points and hi must be inside the unit circle, otherwise the map z{'w) is not 
conformal. On the unit circle we have \dw\ = ^ and 0(w) = ^{\ogz'(w) + log z' (w~^)) , 
where the first and the second term (the Schwarz refiection) are analytic outside and 
inside it, respectively. (Recall that our notation z'{w~^) means dz{u)/du at the point 
u — w~^.) Plugging this into (23), we get: 

■1. , '—1' 
w. 



If ri 

Fi = -TTj— f log^'(w) -dy,\ogz'{w) + 
4:7n J\w\=i 12 



dw ■ 



247ri J\w\= 

1 / 1 -// -udw If , _i.z"{w) 
— — — : f log^^f-u; ) -TTT-. f log ziw )— — -dw 

2Am J\w\=i ^ ^ w 487ri J\w\=i ^ ' z'{w) 

The integrals can be calculated by taking residues either outside or inside the unit circle. 
The poles are at oo, at 0, and at the points and h^. The result is 



Fi = -l|logr^+ E log^-'K"')- E ^ogz'{hf)\ (25) 

^ z'(ai)=0 z>{bi)=<x 



If the potential V{z) is polynomial, V{z) — Yl^=\ tk^^, i-e., tfe = as A; > m for some 
m > 0, then the series for the conformal map z{w) truncates: z{w) — rw -\- uiw~^ 
and 



m— 1 



z'{w) = r — ctiW ^) 

is a polynomial in (all poles 6j of z'{w) merge at the origin). Then the last sum in 
(25) becomes m log r and the formula (25) gives 



m—l 



^i = -^iogk^ n ^ =-^iog|"'n(i-«^«.)l (26) 
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This formula is essentually identical to the genus- 1 correction to the free energy of the 
Hermitian 2-matrix model with a polynomial potential recently computed by Eynard 
[10]. 



4.3 Determinant representation of Fi for polynomial potentials 

For polynomial potentials the genus- 1 correction enjoys an interesting determinant rep- 
resentation. 

Set 

Drr, := det 



\ dtndtidtu I „ , . , , , 

\ ^ J «/ 0<j,fe<m— 1 

Using the residue formula (20) we compute: 

^ _ 1 / dwo f dwm-i det [hj{wj)hk{wj)] , . 

(27ri)"^/|«,o|=i Wo "' J\wm-i\=i Wm-i nS^^'K)^'(wr^) 

Clearly, the determinant in the numerator can be substituted by ^det^{hj{wk)) and 

det [hj(wk)] = (m — 1)! r5"^("*~^)Am(wi), where Am(wi) is the Vandermonde determinant. 
Each integral in (27) is given by the sum of residues at the points inside the unit 
circle (the residues at = vanish). Computing the residues and summing over all 
permutations of the points a^, we get: 

llijjl-'- CljCli) 

As is seen from (10), the last non-zero coefficient of V{z) equals tm — . (We regard 

it as a fixed parameter.) Therefore, HI^i — (—1)'^ m{m — l)r"^''Hrm and we represent 
Fi (26) in the form 

11 1 - 

Fi = — log Dm - — (m^-3m-F3)logr - —(m-1) log const (29) 

where const is a numerical constant. Recalling (18), we see that Fi, for models with 
polynomial potentials of degree m, is expressed through derivatives of Fq: 



1 / d'^F \ 1 d'^F 1 

Fi = — log det ^ ° (m^-3m+3) (m-1) logtm + const (30) 

24 ^rnxm \dtodtjdtk J 24^ ^ dt^ 2A^ J b m V y 



where j, k run from to m — 1 . 

Similar determinant formulas are known for genus- 1 corrections to free energy in 
topological field theories [21]. 



5 Fi from loop equation 

The standard (and powerful) method to obtain 1 /A^-expansions in matrix models is to 
use invariance of the partition function under changes of matrix integration variables. In 
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the 2D Coulomb gas formalism, this reduces to invariance of the partition function (1) 
under diffeomorphisms 

Zi > Zi -\- ai^Zi^ Zi) , Zi > Zi -\- ai^Zi^ Zj^ 
The invariance of the partition function in the first order in e results in the identity 

^ j d,^{e{z,,Zi)e')\{d^z^ = Q (31) 
« i 

for any function e{z^ z). One may read it as Ward idcntitiy obeyed by correlation functions 
of the model. For historical reasons, it is called the loop equation. Since correlation 
functions are variational derivatives of the free energy with respect to the potential, the 
loop equation is an implicit functional relation for the free energy. 



5.1 Loop equation in general normal matrix model 

A closed loop equation does not emerge for the special potential (2). It can be written 
only for the ensemble of normal matrices with a general potential W in (1). Let it be of 
the form 

W{z) = -zz + V{z) + V{z) + U{z) 

where U is only assumed to have a regular Taylor expansion around the origin starting 
from cubic terms. 

Choosing e{zi,Zi) = {z — Zi)~^ and plugging it into (31) with £ given in (5), one is 
able to rewrite (31) as a relation between correlation functions of the field 



$(;2) = -/itr log \{z - M){z- M^) \ = -^^log|^ - 



Zi 



|2 



or d^{z) = —htr{z — M)~^ (here and below d = dz)- Note that d^{z) is trace of the 
resolvent of the matrix M and A$(2;) = —A-np{z), where p is the density of eigenvalues. 
After some simple rearrangings, the loop equation following from (31) acquires the form 

^/^^5^^f^^rf=c = ((»W)^> + ft(a^*W> (32) 

(For any symmetric function /({zj}), the correlation function (/) is defined as the integral 
/ f{{z.i})\/\[^{zi)\^]\jeTi^''^^^(fzj with a normalization factor such that (1) = 1.) This 
relation is exact for any finite N . Supplemented by the relation 

d^{z)) = --^ + d,V{z)F (33) 

(also exact) which directly follows from the definitions of the free energy and the field 
the loop equation allows one to find the ^-expansion of the free energy. 
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5.2 Expanding the loop equation 

The ^-expansion of the free energy for the general normal matrix model is more compli- 
cated than the one discussed in the previous sections. It contains all powers of h, not 
only even: 

log Zjv = Fo + hFif2 + h^F^ + 0{h^) (34) 

so it hardly has a direct topological interpretation. Accordingly, /i-expansions of mean 
values and other correlation functions are expansions in h rather than h^. In particular, 

mz)) = ^o{z) + h^y2{z) + h^^,{z) + 0{h') (35) 

We proceed by expanding the loop equation in powers of h. In the leading order, the 
second term in the r.h.s. vanishes, and z-derivative of the both sides gives: 

{dW{z) - d^o{z)) A^o{z) = (36) 

This just means that for 2; e D the extremum condition (7) is satisfied and A^o{z) — 
otherwise. Inside D, the leading term of the mean density, po{z), is given by po{z) = 
(t{z), where a{z) is defined in (8). (Note that the function cr is defined by this formula 
everywhere in the complex plane, and does not depend on the shape of D, while po 
coincides with cr in D and equals in D"^.) For potentials of the form (2), a{z) = I/tt. 

Being developed into a series in h, the loop equation gives an iterative procedure to 
find the coefficients ^i{z). We need the following results on the correlation functions for 
the general normal matrix ensemble (see [5]): 

d^z)) = / 4^ + 0{h) (37) 
/ Jd C-z 

$(^i)$(^2)) = 2n^ (g{zi, Z2) - G{z,, 00) - G(oo, Z2) - log h^jA\ +o{n') (38) 

where the connected correlation function is defined as (fg) conn ~ (fd) ^ if) id)- Note 
that the function (38) has no singularity at coinciding points ^1,2:2 G D'^. Merging the 
points, we get: 

{id^z)f) ^^{ro;z} + Oih') (39) 

\ / conn (3 

where 

^ w'{z) 2[w'(z)) 

is the Schwarzian derivative of the conformal map w{z). 

After these preparations, further steps are straightforward. Terms of order h and 
of the loop equation give: 

^jL{z,C)A^i/2i0d'C^-d'Mz) 

(40) 



^ I L{zX)A^i{C)d\ = -[{d^i/2{z)y + d'^y2{z) 
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where the kernel of the integral operator in the l.h.s. is 



dW{C,) - d^o{z) 



(41) 



It should be noted that the ^-expansion of the loop equation may break down for 
2; e D. This is mainly because the correlator , when the two points arc close 

to each other and belong to the support of eigenvalues, is not given by cq. (38). At the 
same time, for our purpose we need this correlator just on very small distances, when the 
two points merge. Naively, for z,z' & D the correlator diverges as z' — > z. This means 
that its short-distance behaviour is in fact of a different order in h and must be calculated 
separately. Fortunately, this problem can be avoided by restricting the equations to D*^, 
where no divergency emerges on any scale and one may think that the short-distance 
behaviour of correlation functions is still given by eq. (38). (However, we understand 
that this argument is not rigorous and need to be justified by an actual calculation of 
correlation functions at small scales.) Hereafter, z in (40) is assumed to be outside the 
support of eigenvalues, i.e., the equations should be solved for z e D*^. In this region the 
functions $^(-2) are harmonic. 

From (33) we see that 



The strategy is to find $fc's from (40) and then "to integrate" them to get F^s, i.e., to 
find a functional such that it obeys (42). A general method to find the "derivative" 
T>(z) of any proper functional of the domain D*^ is proposed in [15]. 

An important remark is in order. Suppose we restrict ourselves to the class of models 
with potentials of the form (2) (i.e., with a{z) = l/vr), like in previous sections. Applying 
dzT>{z) to the functional (23), that is Fi in this case, we obtain a wrong answer for 
dz^i{z), which does not obey the loop equation (40)! This seemingly contradicts eqs. (42) 
and so explains why one has to deal with the arbitrary potential. The matter is simply 
that there are functionals such that they vanish for potentials with a{z) = 1/tt but their 
"derivatives", dz'D{z), do not. They do contribute to $1 and restore the right answer. 

5.3 Free energy of the general model 

Skipping further details, we present the results for the general model of normal matrices. 
The answer for Fq is familiar [5]. It is given by (13). The first correction, F1/2, is 



To write down the full answer for Fi in a compact form, we need to introduce, along 
with the 4>{w) (22), another function. 



a,P(z)Fi/2 = a,$i/2(z) , dMz)F, = d,^i{z) 



(42) 




(43) 




(44) 
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and the function x^i^) defined in the domain D*^. It is a harmonic function in D*^ with 
the boundary value The function is the solution to the Dirichlet boundary 

problem: x^{^) — IdD^nG{z,^)x{0\d^\- The explicit formula for Fi reads: 



247r 
1 

^ Stt 



J\w\>l J\w\=l JC 

[ \Vx\'d'z-i xdnX"\dz\-U Axdh 
Jd JdD 2 Jd 



+ 



(45) 



where x in the first three integrals is treated as a function of w through x — xi^iw)). 

The r.h.s. of this formula is naturally decomposed into two parts having completely 
different nature, the "quantum" and "classical" parts: Fi = Fi'^'' + Fi'^^\ The (most 
interesting) quantum part can be again represented through the regularized determinant 
of the Laplace operator in D*^ with Dirichlet boundary conditions. However, now the 
Laplacian should be taken in conformal metric with the conformal factor e^^^^^ . Equiva- 
lently, on the exterior of the unit circle the Laplacian should be taken in the metric with 
the conformal factor e^'^'^'^^^'^^^^'^'^^^ ; we see that and x do enter as the sum + x in the 
first line. More precisely, the formula for regularized determinants of Laplace operators 
in domains with boundary known in the literature (eq. (4.42) in [18]) allows us to identify 

(,) _ 1 det i-e-'^Ac) 

-2^°S det (-e-2xA,.) (^6) 

The classical part comes from "classical" (though of order h) corrections to the shape of 
the support of eigenvalues, which always exist unless a{z) is a constant (see below). It 
is essentially given by F^'^^^ = limn^o{^ {{tr x{M)y) ^^^J . 

Different terms of the ?i-expansion gain a clear interpretation in terms of the collective 
field theory of the normal matrix model, in the spirit of [22]. In this context, it is natural 
to start with the general Coulomb gas model with arbitrary /?. The generalized loop 
equation 

i_ J dwioimo) ^ p{io^z)r) + (2 - m{d^^z)) m 

can be understood as the conformal Ward identity for the collective theory. This allows 
one to find the effective action in the form 

So^pJI p{z) log \z - C|p(C) d'zd'C + I W{z)p{z) d 
S,^{l-^)hj p{z)\ogp{z)d' 



(48) 



z 



The second term, 5*2, is a combination of the short range part — |plogp and the entropy 
plogp. (See [23, 24], where similar actions for unitary and Hermitian matrix ensembles 
were discussed.) 
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This action suggests to rearrange the ^-expansion of the free energy (34) and write it 
in the "topological" form F — Z^g>o T^^Fg, where each term has its own expansion 

F9 = F^'^ + E (49) 

n>l 

Here Up = {2 — (5)h is regarded as an independent parameter. The equilibrium density 
of charges, po, is determined by 5S/5p — which leads to the Liouville-like equation 

-^Alog po{z)+(3po{z)^a{z) (50) 

in the bulk. For (5^2 the first term generates corrections to the shape of the support 
of eigenvalues. The classical free energy is Fq = F^^ 

particular we see that F1/2 given in (43) is in fact F(S^^ while the "classical" part Ff^ of 
(45) is fS^\ The "quantum" part is then F^'^'^ = F^°\ 

The collective field theory approach to the normal and complex matrix ensembles will 
be presented elsewhere. 
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